CELL DENDRITE TRANSITION TIN BASE 


ABSTRACT 


The transition from cellular growth dendritic growth was investigated for 
dilute binary alloys zone-refined tin with lead, bismuth, and antimony 
terms the rate growth the temperature gradient the liquid ahead 
the interface and the solute concentration Co. was found that 
previous relationship describing the transition for lead base alloys applied, 
least first approximation, for the low solute concentration tin results, 
the high concentrations solute, however, was necessary introduce another 
variable, namely the cell size transition d;. The transition relationships only 
applied for alloy systems with distribution coefficient ko, less than unity. For 
systems with (antimony tin) was used which was equivalent 
the method described for determining this ‘‘equivalent 


INTRODUCTION 


Rutter and Chalmers (1953) have shown that when alloy solidifies under 
certain growth conditions, the liquid adjacent the advancing interface 
constitutionally supercooled. The degree constitutional supercooling, shown 
theoretically Tiller (1953), dependent upon the following variables: 
the rate growth the temperature gradient the liquid ahead 
liquid interface and the solute concentration Rutter and Chalmers also 
showed that the cellular growth structure results from small amount 
constitutional supercooling, while Winegard and Chalmers (1954) postulated 
that dendritic growth would occur larger amounts supercooling. This 
was verified experimentally Morris al. (1955). 

The transition from cellular dendritic growth mechanism was investi- 
gated Tiller and Rutter (1956) for alloys terms the growth 
variables. was found that the important parameters describing the tran- 
sition were G/R}, and the orientation the crystal. For given concen- 
tration, the transition occurred over range values, depending upon 
the orientation the crystal. The limits the transition range, referred 
and breakdown, were linear when plotted versus 
The exact orientation effect was not determined except that the upper 
breakdown was the dendrite direction. 

Holmes al. (1957) conducted similar experiments with silver zone-refined 
lead. When these results were compared those Tiller and Rutter, was 
found that the two systems could related expressing the concentration 
Co/ko, i.e., transition occurred the same G/R? ratio for the same Co/ko 
concentration. Tiller (1956) suggested similar relationship and expressed 
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where was constant dependent upon the orientation. The constant 
was postulated function cell size since cell size was dependent upon 


the orientation. 


(2) Ag A'dg 

where the cell size for given orientation, then the general relationship 
becomes 

(3) 


This relationship was not verified experimentally, since cell data was not 


available for lead base alloys. 

The primary purpose this investigation was determine the more 
general equation (equation (1)) applied another alloy system with solvent 
material entirely different crystal structure than lead, namely tin, which 


interlacing body center tetragonal. 


EXPERIMENTAL PROCEDURE 
Three tin alloys were chosen, lead tin, bismuth tin, and antimony 
tin; the distribution coefficients for each alloy are listed Table Antimony 
particular interest since has distribution coefficient greater than 


unity. 
TABLE 
Alloying element Distribution coefficient, 
tin 
0.1 
0.2 


The alloys were prepared from zone-refined tin and chemically pure lead, 
bismuth, antimony. The melting the specimens prior growth was 
done vacuum, but growth was carried out under atmosphere argon. 
All the alloys were not analyzed, since was shown from preliminary experi- 
ments that the percentage weight added agreed with the analyzed value. 

The solidification structures were observed interfaces. The 
interfaces were obtained the use modified Elbaum and Chalmers (1955) 
apparatus, which essentially mechanism for decanting away the liquid 
from the solid during solidification. The apparatus, which two-part graphite 
boat larger graphite container, shown Fig. decanting, 
and after decanting. Solidification began from the water-cooled copper 
block placed the end the movable portion the graphite boat; this 
shown the right side Fig. Freezing was allowed proceed from 
right left until the interface reached the stationary half the boat. The 
movable half was then accelerated the right along the track the larger 
graphite container; the solid was thus pulled away from the liquid expose 


the interface. 
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Fic. view the graphite boat for growing the specimen and for decanting away 
the liquid before decanting, (b) after decanting. 


important when using this technique maintain unidirectional heat 
flow, and two special features were used ensure that this would the 
case. Firstly, the two-part inner boat was surrounded graphite case, and 
the air gap between them helped prevent lateral heat loss. Secondly, 
small graphite insert, seen Fig. between the movable and stationary 
halves the boat, prevented heat being transferred along the graphite boat 
faster than through the solidifying metal. 

The growth conditions were measured manner similar that used 
Tiller and Rutter, except that only two thermocouples were used, and these 
were placed directly the melt apart. The position the thermo- 
couples, which were inside the horizontal alundum tubes, shown Fig. 
The liquid was decanted when the solid had advanced the first thermo- 
couple. Since the thermocouples were connected differential, reading 
this point gave the temperature gradient the liquid directly. 

The growth conditions were varied changing the longitudinal tempera- 
ture gradient tube furnace that surrounded the boat, and varying 
the rate which this gradient was moved along the specimen. This provided 
range values from 1.5 13.0°C/cm and values from 0.001 
0.015 cm/sec; this range larger than that obtained using the thermal valve 


technique Tiller and Rutter. 


| 
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OBSERVATIONS AND RESULTS 


Interface Structures 

The stages the transition from cellular interface dendritic interface 
are shown series photomicrographs Fig. This series was taken 
from 0.2 at.% lead tin alloy, but typical for all concentrations and for 
the other tin alloys investigated. Figure 2(a) the interface 
specimen grown under the highest ratio. Each successive photomicro- 
graph represents interface specimen grown lower G/R? ratios. The 
transition from cells dendrites was defined that point when the cell 
walls became branched disappeared. This shown Fig. 2(d) for this 
particular series. This criterion similar that used Tiller and Rutter, 
and Holmes al. 


Growth Conditions for Transition 

The growth conditions for the transition are shown graphically for the three 
tin alloys Figs. and where the results are plotted Co/ko versus 
The value used for the antimony tin system not but 0.6, 
which defined ko’, the equivalent value. The reason for using this 
value and the method for determining the equivalent discussed later. 
The curves shown the above figures represent the growth conditions for 
transition; the right the curves are the growth conditions for cells, and 
the left are those for dendrites. The transition curves are not linear, and 
there apparent range ratios over which transition occurs for 
given solute concentration. 


Orientation Crystals Observed the Interfaces 

Each specimen examined was polycrystalline with the crystal size dependent 
upon the concentration. low concentrations, the crystal size was very 
large and the orientation each individual crystal could determined. 
random sample the crystals showed that the orientations were within 20° 
the reported dendrite direction [110]. the higher concentrations, the 
crystal size was too small determine the orientation each crystal and 
too large for conclusive determination the preferred orientation. Sufficient 
evidence was obtained, however, indicate that some type preferred 
orientation existed. The structures observed polycrystalline interface 
were similar all the crystals, thus indicating that the solidification structures 
tin were not sensitive small changes orientation. 


Cell Size 

was observed for all three alloy systems that the cell size crystals 
grown the cellular region decreased G/R? decreased. The cell size measure- 
ments, average five readings, are tabulated Table along with the 
corresponding growth conditions. scatter diagram showing the decrease 
cell size decreased for lead tin alloys plotted Fig. the results 
the other alloy systems give similar diagrams. also evident from the 
scatter diagrams that, within the concentration range investigated, the cell 


PLATE 


Transition from cellular dendritic interface for 0.2 at.% alloy 
(150X). Growth condition for each structure is: (a) G/R? 430; (b) 356; G/R? 
200; G/R? 64; (e) G/R? 31. 
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TABLE 
Alloy 
element Co, Cell size 
10.00 285 
9.00 0030 164 00600 
10.75 0021 234 00530 
8.00 0017 195 00560 
7.25 0051 102 00548 
9.50 0005 430 00820 
11.75 356 00553 
11.25 0017 274 00492 
11.00 0017 269 00600 
10.25 0021 223 00525 
7.00 0012 200 00686 
10.25 0028 193 00520 
10.00 0051 141 00522 
7.00 0119 00500 
9.25 0025 185 00615 
8.75 0028 165 00588 
8.75 0051 123 00480 
7.50 0055 101 00520 
8.00 0062 101 00475 
6.50 0068 00425 
6.25 0152 00360 
4.25 0254 00327 
8.50 0038 136 00587 
7.50 0076 00446 
6.00 0152 00408 
4.25 0237 00304 
1.00 10.00 0008 357 00558 
8.00 0019 281 00580 
9.75 0013 271 00604 
9.25 0017 226 00610 
8.75 0019 199 00521 
9.00 0024 184 00615 
10.00 0050 455 00827 
10.75 336 00535 
8.75 243 00617 
208 00497 
9.50 0028 179 00545 
9.00 0042 138 00485 
5.75 
9.00 0021 184 00493 
10.50 152 
7.50 
6.25 
7.50 
5.25 
3.25 
2.4 12.50 391 
1.8 11.00 314 
10.50 284 
9.25 226 
9.75 163 
9.25 149 
7.50 
7.00 
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size independent the solute concentration. The best curves drawn through 
the points for each alloy are shown Fig. these indicate that the cell size 
also independent the solute addition. 

Since the G/R? ratio for transition dependent upon the concentration (see 
Figs. and 5), would expected therefore that the cell size the 
transition point, defined d,, should vary with concentration. This variation, 
shown Fig. for the three tin alloys, linear. order keep the 
curves the same scale, the concentration was expressed Co/ko and the 
value used for antimony tin was again the equivalent 0.6. 
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Fic. Scatter diagram cell size versus the parameter Sn. 
Fic. size versus for all three tin alloys. 
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Cell size Transition, (cm) 
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Fic. Cell size the cell dendrite transition for all three alloys. The equivalent 
value used for the alloy. 


Summary Observations and Results 

(1) The transition from cells dendrites when plotted versus 
was not linear. 

(2) Transition did not occur, for given concentration, over range 
ratios. 
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(3) For cellular growth, the cell size was independent solute concentration 
and solute type, but dependent only upon and 

(4) The cell size transition from cells dendrites varied linearly with 
the solute concentration. 

DISCUSSION 

Comparison Lead and Tin Base Alloys 

The results the lead base alloys and the tin base alloys are compared 
Fig. The lead base alloys are taken directly from the paper Holmes 
al., and the main differences between the two systems are: 

(a) The transition was not found occur over range ratios for 
given concentration. From the crystal orientation measurements, appears 
that least within 20° the dendrite direction there orientation 
effect. 
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5.0 Lead Alloys 


(using 
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Fic. Comparison the cell dendrite transition curves for the tin and lead alloys. 
(Lead alloy curves taken from Holmes al. (1957).) 


The relationship between and for transition not linear for 
the tin alloys; this contrast the results for the lead alloys. The difference 
may part due the fact that the present investigation was carried 
higher ratios. Holmes al., when comparing the two lead alloys, extra- 
polated value while the actual experimental values only went 
3.5 for tin lead and 2.4 for silver lead. The tin results may also 
considered being linear taken only the same value. 
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Since the curves, shown Fig. are close together, seems that 
there major effect solvent crystallography and that the onset 
dendritic growth can predicted from equation (1) least for low con- 
centrations. When the entire concentration range included, however, the 
equation does not apply. This suggests that one the parameters not 
correct that another variable influences the transition. 


Cell Size Transition 
previously shown, the cell size transition decreases the concen- 


tration decreases. This variation can explained terms the cell profile 
and provides further insight into the mechanism the cell dendrite tran- 
sition. Consider the schematic diagram Fig. 10. The diagram represents 
pairs alloys grown various values. The concentration within each 


C,< C2 


for each pair 


(c) 


(d) 


Fic. 10. Effect varying the value the cell profile for two different concentrations. 


pair and respectively, where and the for both con- 
centrations pair the same. The cell size each pair the same since 
was observed that, within the cellular range, cell size was independent 
composition and dependent only upon The diagram shows this decrease 

lowered, may assumed that the amount constitutional 
supercooling ahead the interface tends increase. The increase super- 
cooling can reduced solute taken away from the tip the advancing 
interface and placed the cell wall. This will lead increase the cell 
wall area decrease cell size) decreases. The alloy concentration 
however, has more solute accommodate, and yet has the same cell size 
the alloy composition The extra solute can accounted for the 
cusps the high concentration alloy are deeper than the low concentration 
alloy. Again, this shown Fig. 10. Depending the depth the cusps 


NARA 
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and the cell diameter, the cell profile will tend towards triangular shape 
two dimensions pyramidal shape three dimensions, pyramidal 
instead conical due the platelet system (Rosenberg and Tiller 1957). Once 
the cell profile reaches the pyramidal shape, further decrease cell size does 
not increase the total volume available sink’’. Solute will then 
build ahead the cellular interface, causing larger amounts consti- 
tutional supercooling and thus dendritic growth. This occurs Fig. 10(c) for 
the concentration alloy; any further decrease the cell diameter does not 
increase the total area between the cells for this concentration, i.e., there 
increase area between the cells shown Fig. and (c). This does 
not occur for alloy until lower G/R? ratio reached (Fig. 10(c)). other 
words, for the case, the cells can continue decrease and produce more 
“solute accept the added solute caused the lowering 
since they have not reached the limiting pyramidal shape. 

The above argument not only helps the understanding the cells 
dendrite transition, but also suggests that any relationship between Co, ko, and 
for the breakdown should include the parameter d,, the cell size 


transition. 
(4) 


(where the cell size transition) was applied the three tin alloys 
using the equivalent value for the alloy. 

The three curves are plotted Fig. 11. linear plot was obtained indicating 
that the cell size important variable describing the transition. The 
above relationship cannot tested for the lead base alloys since cell size 
data transition are available. 


Dendrites 
Cells 
o Pb in Sn | 
100 200 300 400 500 


Fic. Growth conditions for the cell dendrite transition with the variable cell 
diameter included. 
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The relationship expressed equation (4) similar that developed 
Tiller (1956) account for the orientation effect observed lead alloys. From 
the present work, appears that cell size even more important than pre- 
viously thought determining when transition occurs for any given alloy, 
since cell size transition depends upon concentration. 


Distribution Coefficient Greater than Unity 
The antimony tin results are plotted Fig. using value and 


equivalent value 0.6. The equivalent value defined the distri- 
bution coefficient type phase diagram which will give exactly the 
same amount constitutional supercooling under the same growth con- 
ditions. find this equivalent first necessary reverse the sign 
the slope the liquidus curve; next, the value solute depletion ahead 
the type alloy must equated against the solute build-up for equiva- 
lent alloy that the values are the same but opposite sign. ko’ 


represents the equivalent ko, then 


The expression used for the solute build-up and depletion are those developed 
Tiller al. (1953) from the work the redistribution solute ahead 
plane interface, assuming the same non-equilibrium conditions. 

The equivalent value for antimony tin 0.6. When the actual 
value used, shown Fig. the dotted line, the results are not 
agreement but they are when the equivalent value used. This because 
the equivalent must used whenever the value used 
indication the freezing point depression such used the constitutional 
supercooling argument. The use equivalent makes all concentration 
readings positive with respect the original concentration. 


CONCLUSIONS 


can concluded that the cell dendrite transition the tin base alloys 
dependent, was found for the lead base alloys, the growth variables. 
This again indicates that the important concept this transition constitu- 
tional supercooling. can also concluded that the general relationship 
describing the transition not markedly affected the solvent crystallo- 
graphy. Although the relationship describing this transition still does not have 
complete theoretical basis, this investigation, however, does point out the 
importance cell size the transition and the effect distribution co- 


efficient greater than unity. 
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THEORY IONOSPHERIC CURRENTS ASSOCIATED 
WITH AURORAE. 


WEAVER? AND SKINNER 


ABSTRACT 


The current distributions induced dynamo action sheet partially 
ionized gas moving with uniform velocity across perpendicular magnetic field 
are derived for the case when the gas contains elliptical region greater ioniza- 
tion density. Special attention paid processes occurring away from the ends 
long narrow ellipse which the ionization density much greater than the 
surrounding region. 

The results are discussed relationship conditions pertaining the 
ionosphere when long homogeneous present. concluded that 
line current, whose magnitude linearly dependent the height-integrated 
Cowling conductivity within the arc, flows along the arc. 


INTRODUCTION 


One the first observed characteristics auroral displays 
close relationship with geomagnetic storms. Although extensive analyses 
geomagnetic disturbances have been conducted statistical and world- 
wide basis, surprisingly little attention has been paid the magnetic distur- 
bances associated with individual auroral displays. 

recent years, has been realized that Chapman’s current systems 
deduced from statistical studies, while providing much valuable information, 
not describe the actual ionospheric currents flowing any particular time. 
particular, Meek (1955) found from his study ionospheric, auroral, and 
magnetic variations that Chapman’s average current systems were inconsistent 
with the individual magnetic disturbances occurring Saskatoon any 
particular time. Therefore, seems preferable examine individual magnetic 
disturbances rather than long-term averages order obtain clearer 
understanding the cause magnetic variations and, particular, their 
relationship with auroral displays. 

this paper, examine the problem from theoretical point view. 
shall accept the hypothesis that dynamo action ionospheric region 
increased ionization density (assumed coincident with auroral display) 
the physical mechanism responsible for induced ionospheric currents, and 
hence those magnetic disturbances associated with the auroral display. With 
this basis, the subsequent work devoted mathematical investigation 
the nature the current systems associated with certain types aurorae. 
The current system thin sheet partially ionized gas containing 
elliptical region greater ionization density calculated. The solutions 
this problem are then adapted for application the ionosphere that the 
current flows associated with long homogeneous auroral arc can found. 
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THE CONDUCTIVITIES PARTIALLY IONIZED GAS 


The conductivity ionized gas the presence magnetic field 
anisotropic. the direction the magnetic field there longitudinal 
conductivity the direction the electric field Pedersen conductivity 
and perpendicular both fields Hall conductivity o2. Thus the current 
density vector can written 


(1) j +o2h XE’ 


where unit vector the direction the magnetic field. and 
denote the respective components the electric field the direction of, 
and perpendicular the magnetic field. 

The anisotropy the conductivity the ionosphere, which can con- 
sidered lightly ionized gas, was pointed out Pedersen (1927) and Cowling 
(1933), who calculated expressions for and respectively. detailed 
examination the possible effect the different conductivities induced 
dynamo currents the ionosphere was undertaken Martyn (1948) and 
Baker and Martyn (1953) attempt account for observed discrepancies 
between the magnitude the field predicted the dynamo theory and 
the observed magnetic variations. 

The conductivities for lightly ionized gas containing electrons, 
positive ions (carrying positive electronic charge), and neutral particles can 
derived most simply using mean-free-path theory. Cowling (1945) showed 
that for the conditions found the ionosphere this gives approximately the 
same expressions the accurate theory. assumed that the density 
electrons equals the density the positive ions, but small compared with 
the density the neutral particles, and also that the positive ions and neutral 
particles have equal masses. can then shown that (Dungey 1958) 


where the electronic charge, and are respectively the masses 
electron and positive ion, and and are their respective collisional 
frequencies with neutral particles. the usual notation and represent 
the gyrofrequencies the electrons and ions respectively defined the 
relation 


has been found that the combination plays important 
role ionospheric electrodynamic equations. Chapman calls this expression 
the Cowling conductivity. Denoting have 
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CURRENT FLOW SHEET LIGHTLY IONIZED GAS 


shall calculate the current system thin sheet lightly ionized gas 
containing elliptical, anomalous patch greater ionization density. The 
gas considered moving with uniform velocity perpendicular 
constant magnetic field. such the problem somewhat artificial nature, 
the physical conditions being idealized for mathematical convenience. 

Consider horizontal plane sheet the partially ionized gas having per- 
fectly insulating regions above and below it. assume that all vertical 
currents are prevented the sheet vertical polarization field set 
the accumulation charge the upper and lower surfaces the sheet. Thus, 
right-handed set axes are defined such that the and axes are the 
plane the sheet and the z-axis vertically upwards the current density 
vector has components j,, and the equation continuity can 


written the form 


provided that sources sinks charges exist the sheet. This equation 
will satisfied introduce current function J(x,y) defined the 


relations 


Let the magnetic field directed vertically downwards perpendicular 
the sheet. Thus (1) —k, where unit vector the direction and 


(1) reduces 


(10) 


where 
(11) 

matrix quantity representing the generalized anisotropic conductivity. 

the ionized gas has uniform mass-velocity there will induced 
dynamo field vXH experienced the gas. If, furthermore, there 
ionization anomaly the gas, charge accumulations will set the 
boundaries due the fact that currents can flow more readily within the 
anomalous region high conductivity than they can the surrounding 
medium. The charge accumulation will increase until has created electro- 
static polarization field which just sufficient retard the current flow within 
the anomaly such extent that the rate arrival charge the boundary 
balanced the rate which charge can leak away the surrounding gas. 
Thus the total electric field can written 


where the potential the polarization field. 
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After substituting for (10) and introducing defined (8), (9) 
have 


(13) 
where 

(15) 
Elimination from these equations gives 

(16) 


where the Laplacian operator two dimensions. Similarly, elimination 
yields 


(17) 


except points the boundary the anomaly, where the net charge density 
longer zero. Hence (15), have also 


(18) 


except the boundary. 

The solutions (18) will derived for the case when the gas contains 
elliptical region greater ionization density. Let the potential the polariza- 
tion field set the accumulation charge the elliptical boundary 
denoted inside the ellipse and outside. Clearly and must satisfy 
the following conditions: 


(i) within the ellipse, 
(ii) outside the ellipse, 
(iii) the boundary the ellipse, 
(iv) finite inside the ellipse, 
(v) the distance from the ellipse approaches infinity, 
(vi) the normal component continuous the boundary the ellipse. 


Let the equation the ellipse 


(19) 
the complex where 
(20) 


Consider the conformal transformation the given 
(21) arc cosh (t/p) 
where 


(22) 
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the interior the ellipse with cut taken along the real axis from 
(24) arc tanh 


The elliptical contour transformed into the line while the remainder 
the contour representing the real axis cut the ¢-plane transformed into 
illustrates the transformation. 

Solutions for and satisfying the conditions (i) (v) are 


(26) 


& 
D 
bid 
Ris 


plane 


Fic. The transformation arc cosh (t/p). 


(23) 
' 
' 
' 
plane 
' 
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where A(w) and B(w) are constants. and can evaluated 
means the boundary condition (vi). This requires that continuous 
Now the ¢-plane 


thus require 


where dashed quantities refer the region outside the ellipse. Substitution 
the partial derivatives and this equation yields, after some lengthy 
algebraic manipulation, the following expressions for and 


when 

(29) A=B=0; 

when 

where 


The angle introduced for convenience the above equation the inclination 

Thus the general solutions for and are 


Transformed into the ¢-plane these are 


(36) 
(37) 


There one-one correspondence between the w-plane and the 
However, two-valued function since plus minus sign could placed 
front the square root sign. one-one correspondence between and 
can set replacing the Riemann surface two sheets with 
the connecting cut from 

special case can let that the ellipse approaches circle. 
When and the axes are rotated that the direction the 
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(38) 
(39) 


(41) 


for and polar co-ordinates: 


where 


and 


These expressions agree with those derived Martyn (1953), who discussed 
the induced current systems and hydrodynamic motions cylindrical blob 
enhanced ionization the ionosphere. 

very special interest also the case when The ellipse then very 
elongated, and approaches strip infinite length and breadth 
When that terms order and higher may neglected have 


will seen that and are order compared with which 


order unity. Hence small compared with which itself approximately 


equal 
The polarization field and current distribution away from the ends the strip 


5 
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(50) 


Terms order and and higher have been neglected the above 
expressions. The appearance separate equations for the two situations 
modulus the complex number positive above the real axis and 
negative below it. 

taking the real parts the expressions for and the polarization 
potentials can found. order determine current distributions, however, 
are more interested the derivatives and These can obtained 
directly from and for 


OS; 0S, dw, 


Thus, 

Similarly, 

hence, 


The current distribution near the origin can now obtained substitution 
(13) and (14). This gives 


and 


particular interest are the perturbation i’, that is, the additional 
currents flowing due the enhancement the ionization within the strip. 
They are defined the equations 


(60) 
(61) 
where 


(62) 
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the current density the ionized gas when contains anomaly. Thus 


have 
(63) 
(64) 

The equations for and not depend explicitly this expected 


since only valid the neighborhood the origin. 
important result which can seen immediately from the equations 


(63), (64), and (65) 


since the only component which not order (6/a). Thus, the dominant 
perturbation current sheet lightly ionized gas which contains very 
long elliptical anomaly greater ionization density linear current flow 
along the ellipse except near the end points. fact, finally let 
that the ellipse becomes infinitely long strip, 


and 


Since the components the generalized conductivity matrix lightly 
ionized gas are all proportional the ionization density (see equations (2), 


(3), and (4)) can write 


2 n 


(69) 


where the constant ratio the ionization density, inside the anomaly 
that outside, n’. Introducing into (67), the expression for the perturbation 


current along the strip becomes 


where and are the Cowling conductivities the gas inside and outside 
the strip. the ionization density inside the strip very much greater than 
that outside (but not great that the gas ceases become only lightly 
ionized) can write and hence 


(71) le o3E;. 


APPLICATION THE IONOSPHERE 


have calculated the current systems thin sheet partially ionized 
gas intersected perpendicularly constant magnetic field. shall now 
attempt modify this purely mathematical configuration that approxi- 
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mates the actual physical conditions found the ionosphere. should 
then possible associate the currents calculated the preceding section 
with current systems flowing the ionosphere. 

For simplicity the ionosphere will considered neutral, lightly ionized 
gas composed electrons, heavy positive ions, and neutral particles having 
the same mass that the positive ions. region auroral luminosity 
one greater ionization density and will therefore correspond the anomalous 
regions introduced into the mathematical model. necessary that even 
within regions auroral luminosity the ionosphere remains only lightly 
ionized; otherwise, the expressions for the conductivities and are 
longer valid. Experimental evidence (e.g. Currie, Forsyth, and Vawter 1953) 
suggests that the maximum ionization density within aurora about 
The density neutral particles the lower regions the ionosphere 
around Thus even within the most highly ionized regions the 
aurora, the atmosphere only lightly ionized. 

Although the representation the ionosphere partially ionized gas 
appears very reasonable the configuration the mathematical model 
does not directly correspond the ionosphere intersected the geomagnetic 
field. Also, certain assumptions which were made the development the 
mathematical theory will require investigation see thev hold the 
ionosphere. 

Following are listed the features the ionosphere which differ from the 
conditions assumed our model: 

(1) The earth and its atmosphere are not horizontal planes. 

(2) The mass-velocity the ionosphere (ionospheric wind velocity) not 
constant; its magnitude and direction change considerably during the day. 

(3) There probably vertical current incident the ionosphere due 
impinging stream solar corpuscles causing aurorae. 

(4) The ionosphere not has finite thickness over 200 km. 

(5) The geomagnetic field not perpendicular the ionosphere; 
inclined angle given the magnetic dip that point. 

(6) ionization anomaly the ionosphere does not necessarily move 
with the ionospheric wind velocity; aurora observed stationary 
move with velocity relative the wind. 

Each the above items will investigated order establish what 
modifications, any, are required before the mathematical results can 
applied the ionosphere. 

The Curvature the Earth 

this work are only concerned with local current systems and their 
related local magnetic disturbance, distinct from world-wide effects, 
are restricting our observations small area the earth’s surface. this 
case the earth and ionosphere can considered horizontal planes good 
approximation. 


Wind Velocity 
Although the ionospheric wind velocity not constant, the only com- 
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ponents its periodic variation which are important are the diurnal and 
semidiurnal ones. The change wind velocity, which according radio 
measurements (Briggs and Spencer 1954) principally semidiurnal, 
negligible during the period time which transient ionospheric currents 
grow and decay. The wind velocity can thus taken having constant 
value even though both the magnitude and direction will have different 


values different times day. 


Incident Vertical Currents 

the gas sheet the mathematical model was assumed that anomaly 
ionization density merely existed; the cause the enhancement was not 
discussed. the ionosphere, however, one must take into account the cause 
aurorae, for according most auroral theories the ionizing agent solar 
corpuscular stream incident the earth’s atmosphere and such constitutes 
vertical electrical current unless the stream electrically neutral. parti- 
cular, such current will affect equation (7) and will also possess magnetic 
field which might responsible for significant part the magnetic distur- 
bance vector measured the earth’s surface. 

The effect the incident vertical current the validity equation (7) 
collisions. this case, incident particles are driven back the lines force 
after coming rest, thereby making contribution horizontal current 
flow. The three dimensional equation continuity clearly reduces equation 
(7) this situation. 

ascertain whether the magnetic field the vertical current important 
not order magnitude calculation was carried out. The current sheet 
diminishing current density was replaced, for mathematical simplicity, 
infinite sheet constant density, the current system being closed infinity. 
This sheet actually has magnetic field greater than that required but this 
immaterial are interested primarily upper bound for the magnetic 
field associated with the solar corpuscular current sheet. Indeed, the solar 
particles accelerate back the lines force out the earth’s atmosphere 
after coming rest, then the effective vertical current and its magnetic field 
are substantially reduced even further. The constant current density the 
sheet was taken amp/cm. This consistent with auroral arc width 
and the estimates either incident electron flux given Vegard 
(1955) incident proton flux given Chamberlain (1954), whichever 
might the primary ionizing agent. was found that the magnitude the 
horizontal magnetic disturbance the surface the earth the neighbor- 
only the order magnitude the magnetic field and therefore negligible 
compared with the much greater magnetic disturbances associated with 


auroral arcs. 


The Thickness the Ionosphere 
was argued Section that currents sheet gas flowed horizontally, 


all vertical flow being impeded vertical polarization field. fact, with 
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perpendicular magnetic field only transverse electric field vXH was effective 
that vertical flow would have been absent any case. Were there 
ionization anomaly the ionosphere, the current flow the ionosphere could 
obtained considering series horizontal sheets, each uniform 
conductivity, lying one top the other. the assumption that the wind 
velocity independent height, the equations would exactly the same 
those for thin sheet with the conductivities replaced corresponding 
height-integrated conductivities, the integrations being throughout the 
conducting regions the ionosphere. 

When the ionosphere aurorally active and therefore contains ionization 
anomaly, the situation complicated the presence the electrostatic 
polarization field caused the accumulation charge the boundary 
the anomaly. again consider the ionosphere series thin layers, 
one top the other, the varying ionization density (and hence conductivity) 
from layer layer causes different charge densities accumulate the 
boundaries different heights, thereby adding vertical component the 
polarization field—grad Thus, since there vertical conductivity 
the ionosphere, current flows the boundary the arc and the assumption 
that div vanishes each horizontal plane longer valid. 

However, order magnitude calculation the charge distribution 
the sides the arc revealed that the vertical current leakage from any 
horizontal layer was least factor smaller than the steady state 
horizontal current flowing the layer. Vertical flow can therefore con- 
sidered negligible, not only the steady state, but during the build-up process 
also. Thus, each horizontal plane div essentially zero, and currents are 
induced each layer the ionosphere according the equations Section 
and can height-integrated directly. particular, the height-integrated 
perturbation current density along auroral arc given where 


The Inclination the Geomagnetic Field 

When the magnetic field makes angle with the plane the gas sheet, 
simple matter (Chapman 1956) show that the correspondingly 
modified forms equations (13) and (14) are 


where 

(75) cos? sin? 
(76) sin 

(77) 

(78) cos? sin? 


Thus the inclination the magnetic field merely involves the introduction 


4 
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the correction factors auroral latitudes 80° and all regions 
the ionosphere above km, Under these conditions, was 
found that 1.03, 1.01, and 1.03. view the fact that 
with presently available numerical data calculations ionospheric conduc- 
tivities involve errors 10%, clear that latitudes where auroral 
arcs occur, only insignificant errors are introduced taking the geomagnetic 
field perpendicular the ionosphere. 


The Motion Aurora 

Until now has been tacitly assumed that the region greater ionization 
density moves with the local wind velocity. fact the equations developed 
Section refer quantities measured co-ordinate system moving with 
the velocity This does not affect the expressions for the current densities 
which remain unaltered when transformed co-ordinate system fixed with 
respect the earth, but care needed transforming polarization fields and 
velocities. the incoming ionizing agent itself (and hence the aurora) 
drifting with velocity the above requirements will met. general, 
however, the velocity auroral motions differs from that the local wind 
velocity (Kim and Currie 

Let consider the physical situation when the region being ionized 
moving with velocity different from that the local wind. For the sake 
argument, will take both velocities the same direction with the 
auroral velocity the greater the two. Consider the trailing boundary the 
moving strip which contains accumulation of, say, negative charge. When 
the ionizing agent displaced the negative charge will left isolated and 
will decay recombining with the positive ions the new boundary the 
ionized strip immediately front them, leaving net accumulation 
negative charge the new boundary. This process will continue the aurora 
advances that the negative charge concentration effectively carried along 
with the boundary the strip. Similar reasoning shows that positive charge 
will carried along with the leading edge that effect, the polarization 
field moves with the aurora, that is, defined with respect reference 
system moving with the velocity the aurora. Thus, referred these axes 
the total electric field S), where the velocity the aurora. 
This transforms into referred axes moving with the wind 
velocity. Thus all our results remain valid for moving aurora providing only 
that our mathematical idealization the boundary line contour still 
acceptable, and that div still zero throughout the region. 

practice, course, the boundary the anomaly will diffuse, the 
concentration charge being spread over finite distance. therefore 
necessary condition that this distance small comparison with the total 
width the strip. The maximum possible velocity the aurora relative the 
ionized particles Thus, the relaxation time the ionization 
decay, the aurora will have travelled distance relative the particles 
before appreciable recombination will have taken place. therefore, require 
the condition satisfied, where the width the arc. 
Measurements auroral movements Kim and Currie (1958) suggest 


1102 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


that the average velocity auroral arcs about 10° cm/min; iono- 
spheric winds are the order Taking average auroral 
arc width, the above condition requires min. This satisfied, far 
can ascertained, recombination reactions taking place within 
aurora. For instance, the time for 90% recombination for the reaction 
about 0.1 sec (Biondi and Brown 1949). therefore 
concluded that the assumption that the boundary moving aurora sharp, 
not too drastic approximation. 

must shown also that div remains zero when the auroral arc has 
velocity this case, referred axes moving with the wind, the charge 
density becomes time dependent taking the form 


(79) 


where the velocity the aurora relative the wind. Thus, 


relation 


(81) 
have 


Order magnitude calculations indicated that |pv’| was least factor 
smaller than hence (81), and equation (82) reduces 
div required for the mathematical theory remain valid for 
moving aurorae. 


SUMMARY 


has long been assumed that line current flows along auroral arc. The 
evidence supplied magnetic records points very strongly this fact. For 
instance, Stagg and Paton (1939) and Heppner (1954) found that the angle 
elevation auroral arc from any point agrees very well with the elevation 
the perpendicular the magnetic disturbance vector measured the same 
point. more striking indication auroral line current the reversal 
the Z-component magnetic disturbance one passes underneath arc. 
This has also been observed Heppner, well many others. 

From theoretical investigation the current distribution partially 
ionized gas containing ionization anomaly, has been shown this paper 
how the auroral line current could produced dynamo action the 
ionosphere. The theory has also provided expression (equation (72)) for the 
magnitude the height-integrated line current density. 

subsequent paper, order magnitude calculations for the theoretical 
line current will presented and some the results predicted the theory 
will compared with observational evidence. 
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THEORY IONOSPHERIC CURRENTS ASSOCIATED 
WITH AURORAE. 


WEAVER? AND SKINNER 


ABSTRACT 


The average ionospheric current flowing along long homogeneous 
calculated numerically the basis the theory developed Part This 
compared with the auroral current necessary produce the observed magnetic 
disturbance Saskatoon (52.1° N., 106.6° W.) and Flin Flon (54.7° N., 102.0° W.) 
for arc which satisfies the conditions required the theory. good order 
magnitude agreement between the two values obtained. 

The theory also shows that, general, the magnitude the line current 
proportional the intensity the arc, but that when the ionospheric wind 
directed along the arc, the line current vanishes. Thus simple explanation the 
occurrence arcs with associated magnetic disturbance provided. 

relationship between the average magnetic disturbance due arcs and the 
ionospheric wind velocity also established. The hourly variation magnetic 
disturbance during aurorally active hours Saskatoon and Flin Flon shown 
predict ionospheric wind velocity variations which agree well with those 
obtained radio measurements. 


INTRODUCTION 
Part (Weaver and Skinner 1960), theoretical expression for the 
ionospheric current associated with long homogeneous auroral arc was 
developed. was found that under the assumed conditions, linear current 
flow was directed along the arc, its magnitude away from the ends being 
given 


where are the Cowling conductivities within and outside the arc, 
respectively, and the horizontal dynamo electric field the x-direction 
the arc. (Since arcs are generally aligned direction, 
shall henceforth refer the x-direction eastward along the arc. The 
y-direction then northward, perpendicular the Important deductions 
can made from this result, and desirable attempt compare these 
with actual geomagnetic and auroral observations. Some the comparisons 
which have already been made data collected Saskatoon (52.1° 
106.6° W.) and Flin Flon (54.7° N., 102.0° W.) are presented this paper. 


NUMERICAL MAGNITUDE AVERAGE AURORAL CURRENTS 
Equation (1) can written the form 


where and are the height-integrated Cowling 
conductivities, the integrations being carried throughout the conducting 
received May 1960. 
Contribution from the Physics Department, University Saskatchewan, Saskatoon, 
Saskatchewan. 
the Department Mathematics, University Saskatchewan, Saskatoon, 
Saskatchewan. 
Can. J. Phys. Vol. 38 (1960) 
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regions the ionosphere. Thus order calculate average value for 
the numerical magnitude first necessary obtain values for and 
from the best available ionospheric data. 

The most straightforward method calculating these values plot 
curves and against altitude and use numerical integration. The 
variation the Cowling conductivity per pair has been 
calculated several authors (e.g. Baker and Martyn 1953; Chapman 1956); 
was recalculated for the present work making use the ionospheric data 
given Bates and Massey (1951) and with the same assumptions those 
made Baker and Martyn. For details refer the reader the papers 
these authors. 

Before proceeding calculate the actual values and various 
altitudes, necessary know the height distribution the ionization 
densities along vertical lines passing through both aurorally active and 
unexcited regions the ionosphere. shall call these distributions N(h) 

recent years rocket flight measurements have shown that the normal 
electron density variation with altitude the ionosphere much smoother 
curve than had been originally assumed. From daytime rocket measurements 
White Sands, New Mexico, Seddon, Pickar, and Jackson (1954) report 
that the ionosphere remains heavily ionized throughout the 100-200 
region. existence sharply defined layers was found. These results have 
been confirmed more recent measurements (Jackson and Seddon 1958). 
adopt here smoothed-out version the electron-distribution curve 
given the 1954 paper assigning the values 10°, 
1.9X and the heights 80, 90, 100, 160, 200, and 250 
respectively. The variation with height between these levels taken 
uniform; thus minor fluctuations the experimental graph are disregarded. 
The adoption this model for open some criticism since based 
data pertaining daytime ionospheric conditions above New Mexico. 
Nighttime ionospheric densities auroral latitudes could differ from the 
chosen values much factor 10. graph the adopted height 

order determine the variation with altitude, the distribution 
ionization within auroral arc must determined. Owing the uncertainty 
the assumptions required for the derivation theoretical distribution, 
assumed here that ionization density proportional auroral light intensity 
and therefore adopt the observed intensity distribution for arcs giving 
also the ionization density distribution. The standard luminosity curve 
Vegard (1921) was used, more recent measurements luminosity (e.g. 
Harang give essentially the same distribution. The base Vegard’s 
curve taken 100 km, the mean altitude the lower limit auroral 
arcs. the numerical value the maximum density arc known 
(corresponding the peak the luminosity curve) the horizontal axis can 
scaled off accordingly units ionization density. 

Information concerning the maximum ionization density auroral forms 
present uncertain state. appears, however, the evidence 
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Altitude (km) 


130 


120 


Altitude (km) 


108 


Fic. Adopted electron density distribution (N’) the 
Fic. Electron density distribution within auroral 


radar measurements, that densities the order exist aurorae (see 
e.g. Currie, Forsyth, and Vawter (1953)); for the present calculation this 
seems reasonable value adopt. (Forsyth (1959) has pointed out that 
the figure probably refers isolated anomalies high ionization 
existing within aurora, and that the mean maximum ionization density 
may smaller.) 

Figure illustrates the adopted density distribution, namely, Vegard’s 
luminosity curve scaled such that the peak ionization density has value 


240 
220 
200 
160 
120 
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and with the cutoff altitude 100 km. The graph represents 
the distribution only within the aurora; above and below the arc normal 
ionospheric conditions, specified Fig. exist. Thus can write 


where measured km. 
The Cowling conductivity distributions and corresponding densities 


and respectively are shown Fig. can seen from the graph that 
the layer which exceeds its peak value only wide, the 


- 
E 
= 
a 
= 
= 
a 


Cowling Conductivity (emu) 


Fic. Variation with altitude Cowling conductivity for normal and aurorally active 
ionospheres. 


peak occurring altitude about 105 km. The effective layer for the 
Cowling conductivity the ionosphere when auroral arc present thus 
extremely thin, and restricted the lower part the Under 
normal ionospheric conditions the maximum value the Cowling conductivity 
occurs approximately the same altitude but reduced thousand fold 


magnitude. 
approximate numerical integration yields 


These values show that and hence equation (2) becomes approxi- 


mately 


(4) 
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will recalled that —v,H where the northward component 
the ionospheric wind velocity perpendicular the arc and the geomagnetic 
field. Thus the horizontal width the arc, the total perturbation 
current flowing eastwards along the arc given 


(5) 


According Briggs and Spencer (1954) the maximum value 
Cambridge, England, during the night hours about cm/sec and since 
has semidiurnal sinusoidal variation its average value approximately 
two-thirds large. Since Saskatoon approximately the same latitude 
Cambridge, adopt the same value for the present calculation. Taking 
0.6T and using the above value for obtain from (5) 


10° amp. 


follows from the discussion the graphs Fig. that this current 
restricted relatively thin region the ionosphere the lower part 
the that near the base the arc. 


COMPARISON WITH AURORAL AND MAGNETIC RECORDS 


not difficult devise method testing whether the theoretically 
predicted magnitude consistent with data collected two stations 
situated the vicinity auroral arc and which are equipped with magneto- 
meters and all-sky cameras. The exact location the auroral arc can 
determined from the all-sky photographs and the magnitude the auroral 
line current which would produce the observed magnetic disturbances can 
then calculated, taking into consideration the effect induced earth 
currents. the assumption that these currents can replaced 
equivalent line current some depth below the surface the earth, some 
simple geometric reasoning shows that 


where subscripts and refer the two stations, and are the distance 
and elevation the arc respectively, and and are the horizontal 
(northward) and vertical (downward) components the magnetic distur- 
bance vector. 

Unfortunately, order carry out such comparison three very stringent 
conditions must satisfied. These are: 

(1) the arc must long and situated such that the magnetic stations are 


near its center; 
(2) the arc must the only form auroral activity present the time 


question; 
(3) apart from the magnetic disturbance associated with the arc, the day 


must be, the whole, magnetically quiet. 
Clearly (1) must satisfied order that the requirements the mathe- 


matical theory fulfilled. Conditions (2) and (3) are necessary one 


1109 


WEAVER AND SKINNER: IONOSPHERIC CURRENTS. 


know, beyond doubt, that the registered magnetic disturbance associated 
with the observed arc. Thus, although the method simple, the number 
arcs which lend themselves analysis this manner very small. 

Magnetometers and all-sky cameras have been operating Saskatoon and 
Flin Flon for number years. Both stations lie near the auroral zone and 
are only about 400 apart that the same auroral arc often recorded 
both cameras. examination the records from these two stations 
revealed arc for which conditions (1) (3) above were well satisfied. This 
occurred 0931 G.M.T. March 14, 1956. 

Figure illustrates the essential steps the analysis this arc. represents 
map the area with the position the arc (assumed altitude 


49°N 


Fic. Geographical location the arc March 14, 1956. 


100 km) sketched in, and showing also the magnetic disturbance vectors 
registered each station. Substitution the relevant quantities equation 
(6) yielded 

amp. 


This value agrees extremely well with the predicted order magnitude for 
the average auroral current 10° amp. The fact that somewhat smaller 
exactly what one would expect since the influence the return current 
has been completely ignored the theoretical treatment where was assumed 
the arc was infinitely long. For finite length, the return current densities are 
small compared with the auroral current (they are order see equation 
(64), Part but over large distances their integrated effect becomes important. 

not easy make extensive comparisons this type owing the 
stringent conditions required. However, the theory can tested more 
indirect way considering other predictions which can extracted from 


equation (5). 


60°N 
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ARCS WITH MAGNETIC DISTURBANCE 


The occurrence arcs with little associated magnetic disturbance has 
been noticed since interest auroral studies first began. Vegard (1916) drew 
attention this fact over years ago, and observed that even arcs with 
relatively strong luminosity could possess this characteristic. concluded 
that aurorae general give rise only relatively small magnetic effects and 
are not directly related strong disturbances. However, investigations 
Harang Murcray (1958), and others indicate that increase 
auroral luminosity leads increase magnetic disturbance, fact which 
contradicts Vegard’s conclusions. 

The apparent discrepancy can very simply explained. Equation (5) 
shows that, general, the auroral current indeed proportional the 
auroral intensity, since directly proportional the ionization density 
within the arc, and hence its luminosity. However, obtain 
irrespective the magnitude K;. Thus, even very bright arc could 
produce zero magnetic disturbance the ionospheric wind happened 
directed along the arc that time. The appearance such arcs therefore 
not strange phenomenon; natural consequence the dynamo theory 
for auroral currents. 

will noticed that the explanation offered here applies only auroral 
forms whose lengths are large compared with their widths. therefore 
interesting note that Vegard observed only arcs, bands, and draperies 
possessing this feature although generalized his conclusions all types 
aurorae. 

was not possible test this explanation experimentally since information 
(apart from average results) the direction the ionospheric wind the 
time occurrence arc was available. However, examination the 
records showed that auroral forms without associated magnetic disturbance 
were means rare occurrence. Seventeen such cases were found the 
Flin Flon records for the period December 1955 July 1957 and 
Saskatoon for the period December 1955 April 1957. (For the purposes 
this examination any disturbance less than was considered negligible.) 


MEAN MAGNETIC DISTURBANCE 


Consider long homogeneous auroral arc situated distance from some 
station. Let the orientation the arc such that the angle between geomag- 
netic north and the direction the arc measured clockwise direction. 
Then the magnitude the magnetic disturbances vector AF, directed north- 
wards sense, given 


(7) 

where given (5). Thus, the horizontal component this 
disturbance 

(8) 
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where the height the aurora. That is, 
(9) AH(t) 


where and the variable has been introduced indicate which 
quantities are dependent local time. the northward component 


AH(t), have 
(10) AX(t) sin 
now assumed that the variables v,, and are statistically independent. 


Thus taking the average over many arcs occurring different days the 
time and noting that and sin must independent time, (10) can 


written 


(11) AX(t) 


Our assumption statistical independence may not entirely true, for 
although Jensen and Currie (1953) report systematic diurnal variation 
the orientations arcs, there the other hand, reason suppose that the 
latitude arcs correlated with local time (Vegard 1912). 

Equation (11) implies that the average X-component the magnetic 
disturbance field due auroral arcs varies with local time the same way 
the southward component the average ionospheric wind velocity. 

order test this conclusion, analysis the mean magnetic distur- 
bance for those hours the night during which auroral arc was present 
was made. The maximum value the X-component magnetic disturbance 
for each hour auroral activity (in the form long arcs) was recorded for 
the 1956 and 1957 both Saskatoon and Flin Flon stations. The 
magnitude the disturbance was indicated K-index which was devised 
Meek (1956) for Saskatoon; this index also suitable for Flin Flon. Since 
the variation this region only the order its effect was ignored, 
and the deviation the trace the magnetometer chart from its normal 
position was taken direct measurement the disturbance associated 
with the aurora. (The magnetometers actually record the geographic, not 
the geomagnetic components, but the difference slight.) 

The variation with local time the mean K-index (averaged for each hour 
over the 2-year period) shown histogram form Fig. for both Saskatoon 
and Flin Flon. Smooth curves have been superimposed the histograms. 

Equation (11) predicts that this variation should directly proportional 
the variation the southward component the ionospheric wind. The 
idealized semidiurnal variations the north-south ionospheric wind, based 
experimental data obtained radio methods Cambridge (Briggs and 
Spencer 1954) and Ottawa (Chapman 1953) are also shown for comparison 
Fig. 

Although the nighttime parts these curves are essentially extrapolations 
daytime measurements, they are confirmed limited number measure- 
ments ionization movements made during the night. can 


1112 CANADIAN JOURNAL PHYSICS. VOL. 38, 1960 


Saskatoon 


Mean K- index 


' Flin Flon 


Wind velocity 


Idealized variations with local time the mean K-index and north-south component 
the ionospheric wind. Briggs and Spencer (1954); Chapman (1953). 


seen that these curves are consistent with the nighttime wind variations 
predicted equation (11) from the variation the mean magnetic index. 
The average wind velocity southward before midnight when the mean 
magnetic disturbance positive, and northward when negative after 
midnight. not possible, course, compare the amplitude the curves 
without knowledge the distribution the geographical position the 
arcs, but has already been established Section that there good order 
magnitude agreement, least for the particular arc examined. 

The slight phase difference about hours between the curves not 
considered serious the radio measurements ionospheric movements during 
the night are not sufficiently numerous establish the precise hour reversal 
the wind. 


SUMMARY 


this paper, further discussion the dynamo theory for auroral currents 
(developed Part has shown that the theory explains several observed 
auroral phenomena. particular, the following results have been noted. 

The theory explains the line current flowing along auroral arc. 

Calculations have shown that this current generally confined within 
region less than vertical extent from the base the arc. 

general, the intensity magnetic disturbance associated with arc 
proportional the intensity the arc. 
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times when the ionospheric wind directed along the arc, however, 
there will associated magnetic effects, regardless the intensity the 

The average horizontal north-component magnetic disturbance 
associated with arcs varies time with the south-component the local 
average ionospheric wind. 

Although the conditions involved the mathematical development pre- 
vented extensive comparison these results with recorded magnetic and 
auroral data, for the one particular arc examined there was good order 
magnitude agreement between the theoretical auroral current and the 
observed associated magnetic disturbance. Also, examination magnetic 
data showed that average wind variations were consistent with those found 
experimentally radio methods. When better experimental knowledge 
the behavior ionospheric winds known, should possible test the 
theory more fully. 
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SCATTERING RADIO WAVES IONIZED GAS 
THERMAL EQUILIBRIUM! 


FEJER 


ABSTRACT 


theory developed for the scattering radio waves density fluctuations 
which exist ionized gas thermal equilibrium. 

Expressions for the frequency power spectrum the scattered waves are 
obtained. These expressions make possible interpret the results observa- 
tions this type scattering from the ionosphere terms electron density 
and temperature. 

shown that the characteristic scale the scattering irregularities (this 
scale depends the wavelength the incident radio wave and the scattering 
angle) much greater than the Debye length then the width the spectrum 
the scattered signal determined the thermal velocities (and the collision 
frequencies the latter are sufficiently high) the positive ions, rather than 


the electrons. 
the characteristic scale greater than the Debye length then for low 


collision frequencies the spectrum flat-topped, with two slightly raised shoulders 
situated symmetrically above and below the frequency the incident wave. For 
high collision frequencies the spectrum has only one maximum situated the 
frequency the incident wave. 


INTRODUCTION 


Gordon (1958) has recently pointed out that scattering radio waves 
from ionized gas thermal equilibrium may detected powerful 
radar. His analysis was based the assumption that the electrons were 
random thermal motion the same type the motion executed the neutral 
particles. this assumption was able calculate the scattering cross 
section and the power spectrum scattering. 

Bowles (1958, 1959) finds that powerful radar beamed vertically up- 
wards, weak noiselike return observed from ionospheric heights. Bowles’ 
observations show that the amount power contained these backscattered 
waves approximate agreement with Gordon’s predictions but that the 
power spread over much narrower frequency band than Gordon predicted. 

Bowles (1959) also gives tentative theoretical explanation his experi- 
mental results. His explanation makes use the work Pines and Bohm 
(1952), who considered the density fluctuations electron gas the 
presence uniform positive charge density. Bowles (1959) points out that 
the positive charge density not uniform but fluctuates result the 
random thermal motion positive ions. The electrostatic forces considered 
Pines and Bohm suppress the fast spontaneous fluctuations electron 
density but the same electrostatic forces cause slow fluctuations electron 
density which neutralize the slow spontaneous fluctuations ion density. 
This explains the observed narrow spectrum. 
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Section this paper the treatment Pines and Bohm modified 
removing the restrictive assumption uniform positive charge density. 
The scattered power then calculated with the aid the revised treatment. 
The results this analysis agree with the conclusions Kahn (1959), who 
used slightly different methods his calculation. 

Sections and the present paper the power spectrum the scattered 
signal calculated for number different approximations. 


CALCULATION THE SCATTERED POWER 


Consider medium which the mean dielectric constant Assume 
further that there are irregular deviations from the value within finite 
volume Then very long distance from the volume the electric 
field the waves scattered these irregularities from incident plane 
wave, the electric field which r)] given Booker and 
Gordon (1950) 


where the circular frequency, the time, the radius vector from the 
origin the co-ordinate system, which taken inside the volume 
dxdydz the volume element, the wave number the medium, 
and are the propagation vectors the incident and scattered waves re- 
spectively (the propagation vectors the waves scattered from different points 
within the volume are almost identical because the scattered field cal- 
culated for point sufficiently far from V), and the angle between the 
electric field vector the incident wave and the direction scattering. The 
rationalized system units used. 

ionized gas sufficiently high frequency 


where the number density, the mass, and the charge electrons. 
The deviation the dielectric constant therefore given 


Equation (1) may then written approximately 


where the velocity electromagnetic waves free space. Equation 
(4) shows that the scattered field proportional that spatial Fourier 
component the electron density deviation, AN, whose vector wave number 
is 

This spatial Fourier component itself undetermined but the mean square 
its absolute value may calculated considerations statistical mechan- 
ics. Such calculation, based the assumed validity the 
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mann energy distribution, carried out the Appendix; closely follows 
the method indicated the Appendix Pines and Bohm (1952), but discreet 
singly charged positive ions are assumed rather than uniform positive charge 
density. The more general case positive ions which carry whole multiple 
the elementary charge treated Kahn (1959) slightly different 
methods. 

Equations (76) and (78) the present Appendix may combined 
give 


where Boltzmann’s constant and the absolute temperature. This 
equation shows that the mean square value the Fourier component 
defined the left side equation (5) equals when 
For the mean square value the Fourier 
component equals This conclusion differs from that Pines and Bohm 
mainly that the Fourier component reduced only factor 
for and does not become vanishingly small. The 
difference due the assumptions made the present analysis about the 
discreet nature the positive charge. 

The combination equation (5) with equation (4) leads 


Since the ratio the power density the scattered wave the 
power density the incident wave, one may write for the scattering cross 
section defined the power scattered per unit solid angle, unit incident 
power density, and unit scattering volume 


where the relationship has been used. more convenient 
write equation (7) terms the Debye length the wave- 
length the radio wave the angle between the propagation vec- 
tors the incident and the scattered waves, and the scattering coefficient 
sin single electron. With the aid the identity 


The characteristic length which occurs equation (8), has been intro- 
duced into the study electrolytes Debye and Hueckel (1923); for short 
discussion its significance plasma physics the reader referred Spitzer 
(1956). 

should noted that the derivation equation (1) for the scattered 
field and equation (8) for the scattering cross section was tacitly assumed 
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that the irregular deviations from the mean value the dielectric constant 
not vary with time. The results remain valid very high degree 
approximation, however, the “apparent the irregularities 
(which the present case are the order the thermal velocities the 
particles) are much smaller than the velocity the radio wave. This condition 
usually satisfied practice and equation (1) used the present paper 
calculate the fluctuations the scattered field caused the fluctuations 


Ae. 
For sufficiently small values equation (8) reduces 


(9) 
possible show simpler and more fundamental methods than those 
the present Appendix that this value the scattering cross section results, 
quite apart from any assumptions the velocity distribution, the spatial 
probability density electron constant throughout the scattering volume 
and the electrons are distributed independently space. Such distribu- 
tion electrons will called the present paper. The above simple 
result for random distribution scatterers has been derived Rayleizh 
(1871) and somewhat greater detail Kahn (1959). The problem essen- 
tially the same the problem shot effect temperature-limited diode 
where the probability density time electron constant and the 
electrons are independently distributed time; the above result could there- 
fore also derived methods used Rice (1944). 

One may then state that for sufficiently small values the radio wave- 
length the electrons scatter they were randomly distributed. 

For sufficiently large values equation (8) reduces Thus 
this case the scattering cross section only half what would the 
electrons were randomly distributed. The transition occurs approximately for 


NATURE THE PROCESSES DETERMINING THE FREQUENCY 
SPECTRUM THE SCATTERED POWER 

the previous section the total scattered power has been calculated but 
information its frequency distribution has been obtained. instruc- 
tive consider the physical processes determining the spectrum the 
scattered power semiquantitative manner before the actual calculation 
carried out the next section. 

ionized gas there generally irregular, fluctuating electrostatic 
field determined the irregular distribution the charge density. Such 
electrostatic field tends redistribute the charged particles 
that they satisfy the energy distribution. the poten- 
tial the electrostatic field the field tries redistribute the charged 
particles with charge that their number density becomes proportional 


(10) 
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for the moment one ignores the field set result redistribution. 
the corresponding irregularities density are large as, larger than the 
original irregularities responsible for the electrostatic field, then likely 
that the effect the electrostatic field must taken into account. 

These considerations may expressed semiquantitative manner 
considering Fourier component the electron density 
fluctuations existing given moment and assuming that the positive ions 
are uniformly distributed. The potential the electrostatic field caused 


(12) 


This potential tries redistribute the electrons, changing their number 
density 


(13) 


then this change, which under the influence 
the field tries establish, reduces substantially the original deviation 
the number density. This the conclusion Pines and Bohm (1952), who 
showed that random irregularities larger dimensions than the Debye length 
tend smoothed out. The present arguments are more qualitative 
nature than those Pines and Bohm, which the Appendix the present 
paper based. 

One may express this result Pines and Bohm somewhat differently 
regarding the actual electron density deviation the sum the spontaneous 
deviation which would have existed the absence the electrostatic field 
and the deviation caused the electrostatic field, ignoring temporarily the 
fact that the electrostatic field actually caused the electron density 
deviations. Thus then the fluctuations electron density 
caused the electrostatic fields almost exactly cancel the spontaneous 
fluctuations. These spontaneous fluctuations would caused the random 
thermal motion the electrons the absence the electrostatic field. 

actual gas the density positive ions not uniform. The fluctuation 
ion density may also regarded the sum the spontaneous deviation 
and the deviation caused the electrostatic field. view the very much 
greater mass the positive ions the rate spontaneous fluctuations ion 
density very much slower than the rate spontaneous fluctuations 
electron density. therefore reasonable assume that actual gas 
the electrostatic field has fast and slowly fluctuating component. The 
fast fluctuating component causes the almost complete cancellation the 
fast spontaneous electron density fluctuations but has only 
effect the motion the positive ions. The slowly fluctuating field affects 
the motion both ions and electrons and ensures that the slow electron 
density fluctuations (caused these fields) almost neutralize the slow ion 


| 
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density fluctuations the results the Appendix show. further reason- 
able assume that the electrons, account their low inertia, have ample 
time assume their equilibrium distribution, indicated equation (10) 
corresponding the instantaneous value the slowly varying component 
the electrostatic field. the instantaneous distribution electron density 
almost identical with the instantaneous distribution ion density, the 
slowly varying component the electrostatic field always assumes con- 
figuration which, given time, would change hypothetical uniform ion 
distribution distribution with exactly opposite ion density deviations 
those actually existing the gas. This explains qualitatively why the electro- 
static fields reduce the mean square value the Fourier component the 
ion density deviation for sufficiently small wave numbers factor 
shown the results the Appendix. 

the relatively small effect the electric field the motion the ions 
were neglected, then the calculation the spectrum the scattered power 
would relatively simple. For sin the spectrum would 
that appropriate electrons executing their random thermal motion the 
temperature the gas. For sin the spectrum would that 
appropriate ions executing their random thermal motion corresponding 
the temperature the gas, provided that these ions are assumed scatter 
radio waves they were electrons. This last statement is, course, not 
quite correct view what has already been said about the effect the 
electrostatic field the motion the ions; this effect examined greater 


detail the next section. 
FREQUENCY SPECTRUM THE SCATTERED POWER 


4.1 Introduction 
The following calculation the frequency spectrum the scattered power 


carried out two steps. first assumed that the scatterers execute their 
random thermal motion uninfluenced fluctuating electrostatic fields. Only 
the two extreme cases, which the mean free path the particles either 
very large very small are considered. The results this calculation are 
strictly valid They are valid for 
sin provided that the mass which determines the thermal 
motion the particles taken the ionic mass but the mass which deter- 
mines the scattering properties taken the electronic mass. 

the second part the calculation this latter result 
corrected for the effect the electrostatic field the motion the positive 
ions. The same two extreme cases are again considered. 


4.2 Scattering Charged Particles Random Thermal Motion 

the effect collisions negligible then the scattering particles may 
sorted into groups such that the velocities all particles belonging the 
same group are almost identical and that the particles within each group are 
randomly distributed space. The power scattered single group is, 
then, terms equation (9) and the subsequent discussion, equal the 


4 
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product the scattering coefficient single particle and the number 
particles the group. The frequency the signal scattered single group 
differs from the frequency the incident wave the Doppler shift caused 
the common velocity the particles. This suggests that more con- 
venient sort the particles into groups such that the velocities particles 
within the same group lead the same Doppler shift. The shape the fre- 
quency spectrum then determined the distribution the particles into 
these groups and thus the velocity distribution the particles. 

The Doppler shift the frequency the scattered wave caused that 
component the particle velocity which perpendicular the bisector 
the angle formed the direction incidence and the direction scattering. 
written for this component, the Doppler shift 


where and are the circular frequencies the incident and the scattered 
wave respectively. 

The Maxwellian distribution velocities, written terms single 
component, shows that the number particles with velocity components 


(15) 


and therefore, the value from equation (14) substituted, the number 
particles which scatter with Doppler shifts between and w+dw—wy 
given 


(16) Si(w) exp[— 


The quantity also proportional the power scattered within the 
frequency range extending from view the introductory remarks 
this section. Therefore, the effect collisions negligible, then the 
frequency power spectrum has the Gaussian form indicated 
equation (16). The function also shown the broken curve Fig. 

Collisions have the effect ‘‘slowing down” the motion individual 
particle continuously changing the direction its travel, and they may 
thereby reduce the spectral spreading just calculated. The fraction particles 
which, time interval long compared the collision interval, suffer dis- 
placements which have components given direction between and x+dx 
proportional (see, for example, equation (172) Chandrasekhar (1943)) 


(18) exp(—x?/4Dt)dx 


where the diffusion coefficient the particle. 
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RELATIVE SPECTRAL POWER DENSITY 


(w 


Fic. The spectrum functions and the low collision frequency approxima- 
tion. given equation (17). 


specifying the particle distribution the beginning and end the time 
interval one free assume that the particles have moved with hypo- 
thetical velocity distribution 
(19) 
and such velocity distribution would the absence collisions lead 
frequency spectrum 

This frequency spectrum does not actually occur, but the relationship between 
the particle distributions existing the beginning and end the time interval 
can determined correctly applying it. The hypothetical autocorrelation 
function obtained taking the Fourier transform the power spectrum 


(20) therefore gives the correct autocorrelation coefficient the received 
signal for This hypothetical autocorrelation function given 


The value applicable the actual problem only for the one value 
but the calculation may repeated for each turn. This then leads 
the correct autocorrelation function the physical problem 


The Fourier transform (22) proportional the frequency power spectrum 
the scattered waves. given 
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The function shown the broken curve Fig. 
The coefficient diffusion may written the form (Chandrasekhar 
1943) 


(24) KT/mB 


where the frictional coefficient the Langevin equation for the motion 
the particle, which for the present purposes may treated Brownian 
particle. the particles were electrons and collisions were mainly with neutral 
particles then would equal the electron-neutral collision frequency. 


° 


RELATIVE SPECTRAL POWER DENSITY 
° 


-4 -3 = 2 3 4 5 


Fic. The spectrum functions S2(w) and So2(w) the high collision frequency approxima- 
tion. given equation (38). 


The spectrum functions given equations (16) and (23) represent two 
extreme cases. Some rough criterion needed when collisions are 
sufficiently frequent modify equation (16) and make equation (23) better 
approximation. Collisions are expected make the spectrum narrower and 
therefore would reasonable assume that (23) better approximation 
than (16) when the power band width the spectrum expressed 
narrower. 

Integration shows that the band the spectra expressed 
(16) and (23) are equal for 


where 


the mean free path the particle. 
one would therefore use equation (16) for the spectrum; 
equation (23) would appropriate. mentioned the 
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introduction the present section, the results this subsection are strictly 
and (23) must taken the values appropriate electrons. 

the results the present subsection may still regarded 
first approximation the power spectrum scattering provided that 
and represent values appropriate ions. Thus the first approxi- 
mation the width the spectrum very considerably reduced when 
but the shape the spectrum not changed very much. 
The more correct treatment the next subsection shows that, the effect 
the electrostatic fields the ion motion taken into account, this results 
considerable change the shape the spectrum addition the reduced 
width when sin 


4.3 Effect the Electrostatic Field the Spectrum 

this subsection the condition assumed and the approxi- 
mate results the previous subsection are corrected for the influence the 
electrostatic field, caused the momentary charge distribution, the 
motion the ions. 

The high collision frequency approximation, characterized sin 
with /,, appropriate ions, considered first. The first 
the spectrum scattered power then given S2(w) equations (23) and 
(24) with values and appropriate ions. Calculation the correction 
the spectrum function for the effect the electrostatic field the 
motion the ions follows. 

was shown Section that the amplitude the wave scattered 
given direction determined the spatial Fourier component with wave 
number the electron density deviation AN. The variation with time 
the scattered field may therefore determined from the variation with 
time the Fourier component, which according equation (4) proportional 
The spatial Fourier component may written the form 
Fourier integral time (omitting numerical factors involving the 


form 
(27) 

Since proportional this expression one has 
(28) exp 


Equation (28) shows that proportional the amplitude 
those scattered waves whose frequency whose frequency 
differs w—wo from the frequency the incident wave. 

electron density deviation. The total electron density deviation may 
written terms the Fourier component (27) the form 


; 
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that elementary contribution this integral the vicinity 
and may written the form 


The effect the electrostatic field the ion motion causes the actual ion 
number density deviation differ from the hypothetical ion number 
density deviation 6V, which would have existed the spontaneous random 
thermal motion the ions had not been modified the electrostatic field. 
One may regard the sum this spontaneous value and the deviation 
caused the electrostatic field. This electrostatic field, the slowly 
varying component Section has just the correct value keep the electron 
density deviation equal the ion density deviation. Since the ion density 
corresponding the deviation the electrostatic potential 
required make the electron density also equal satisfies relation 
which follows from equation (10). This relation 


(31) 
(32) 

The deviation caused this potential may determined from the 
diffusion equation (see, for example, Chandrasekhar (1943) equation (311)) 
(33) (Ne/Bm) div grad div grad 


and with the aid equations (30) and (32) one may rewrite equation (33) 
the form 


and thus 


Equation (35) shows that the Fourier component the spontaneous electron 
density fluctuations wave number reduced factor 


and thus the power density the scattered waves the vicinity the 
frequencies and reduced factor equal the square 
expression (36). Using equation (24) and the identity 
the reduction factor may written the form 

Abbreviating 


the spectrum described (23) takes the form 
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(39) 
and after multiplication the correction factor (37) the corrected spectrum 
function becomes 
(40) 
The function shown the solid curve Fig. Integration shows 


that 


agreement with the theory Section The reduction the total scattered 
power factor caused mainly the reduction the power scattered 
near the frequency the incident wave, the solid curve Fig. shows. 
Equation (40) applies the high collision frequency approximation for 
where has the value given equation (26) for positive 


For low collision frequencies sin with appropriate ions) 


the first approximation the spectrum the scattered power given 
equations (16) and (17) with taken the mass ion. Calculation the 
correction these equations closely follows the calculation for high collision 
frequencies equation (32). The deviation caused the electrostatic 
field not determined, however, the diffusion equation. The appropriate 
equation which governs the motion positive ions the electrostatic field 
derived from the potential equation (249) Chandrasekhar (1943) with 


the right-hand side set equal zero 


(42) 
where the distribution function positive ions (over positions and 
velocities), grad, calculated differentiation with respect spatial 
co-ordinates, and grad, calculated differentiation with respect 
velocity components. the force per unit mass acting the ions and the 


present case given 


the velocity distribution the unperturbed gas assumed Maxwellian, 


one may write 


where the perturbation function the velocity. this expression 
substituted into equation (42), one obtains for the component after 


linearization 


(45) 
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Substitution from (43) into (45) leads 


and with the aid (30) one may rewrite equation (46) the form 


Uz 
|k—kojv, 


One has integrate the perturbation the 
distribution function over all velocity space obtain the fractional per- 
turbation the ion number density caused the electrostatic 
potential Only the integration over the component has carried 
out detail, the other two components merely leading toa factor 
With the aid equation (47) one obtains for the perturbation ion density 


(47) 


The integrand the right equation (48) has infinity 

The evaluation this type integral has been first discussed connection 
with plasma oscillation Landau (1946). The results are reviewed Berz 
(1956). terms these results the integral equation (48) can replaced 
its principal value plus times the residue the pole the integrand 
(cf. equation Berz). This leads 


0 


where defined (17) with denoting the mass positive ion. 
Equation (49) may rewritten with the help the relation 
the form 


ww’ 
0 


Equation (50) shows how the spontaneous fluctuations 6.V, are modified 
the action electrostatic fields the absence collisions. Since the scattered 
power proportional the square the electron density deviation and since 
the uncorrected spectrum given equation (16) the corresponding 
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This spectrum shown the solid curve Fig. Numerical integration 


shows that 


agreement with the theory Section 2.* The reduction the total scattered 
power factor caused again mainly reduction the power 
scattered near the frequency the incident wave. The resulting spectrum 
function flat-topped with slightly raised shoulders symmetrically placed 


about the frequency the incident wave. 


CONCLUSIONS 


expression for the scattering cross section has been obtained the 
form equation (8) terms the scattering coefficient single electron, 
the number density electrons, the wavelength the incident wave, 
the angle between the propagation vectors the incident and scattered 
waves, and the Debye length then equation (8) results 
the scattering cross section No, while sin(@/2) equation (8) 
yields 

short review the results obtained for the frequency spectrum the 
scattered power let the case considered first. the 
same time where the mean free path electrons, the 
spectrum has Gaussian distribution given equation (16) and (17) with 
indicating the mass electron. the other hand then 
the spectrum narrower and given equation (23) with taken the 
diffusion coefficient electrons. 

Considering the case one finds that the spectrum very 
much narrower because its width largely determined the thermal motion 
the ions rather than the electrons. For sin(@/2) the spectrum 
function given approximately equations (40) and (38) with taken 
the diffusion coefficient ions. The shape the spectrum function shown 
the solid curve Fig. 

For sin(@/2) the spectrum function given equations (51) and 
(17), with taken the mass ion, obtained. The shape the 
spectrum function shown the solid curve Fig. 

The effects uniform magnetic field the scattering radio waves 
ionized gas thermal equilibrium were not considered the present paper. 

*Equation (51) has been derived independently Dougherty and Farley 
the Cavendish Laboratory, using different methods; they also establish equation (52) analy- 
tically (private communication). 
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is, however, reasonable assume that if, certain ionospheric applica- 
tions, the gyroradius positive ions much greater than and 
not nearly perpendicular the magnetic field then the results 
the present paper are not seriously affected the presence the magnetic 
field. If, however, these conditions are not fulfilled, reformulation the 
theory becomes necessary. 
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APPENDIX 
FOURIER ANALYSIS THE SPONTANEOUS DENSITY FLUCTUATIONS BASED 
THE ENERGY DISTRIBUTION 


The present calculation closely follows that given Appendix paper 
Pines and Bohm (1952), who obtain the mean square amplitude the 
spatial Fourier component electron density fluctuations. Pines and Bohm 
assume background uniform positive charge density while the present 
calculation assumes discreet singly charged positive ions. 

Attention should drawn this stage apparent error the calcula- 
tion Pines and Bohm (1952). their equation (A6) the second term inside 
and therefore the expression the right the final equation should 
and not 

gas containing the average electrons and singly charged positive 
ions per unit volume considered. unspecified number neutral particles 
may also present. Thermal equilibrium assumed. further assumed 
that the probability that the particle lies between and (i.e. that 
its three Cartesian co-ordinates lie between and and 
and while its momentum lies between and (the notation 
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where dr, denotes the volume element and the total energy the 
system, the sum kinetic and potential energy, given 


n m+n 


Here e,, and are the charges the mth and charged particles, the 
distance between them, the mass the nth particle, and the dielectric 
constant free space. the integration over momentum space carried out, 
the probability given spatial distribution the charged particles 
obtained the form 


n m+n 


Let the available space divided into cells equal volume such that 
each cell contains large number charged particles and yet each cell 
small enough that important physical property changes within cell. 
Let the 7th cell contain singly charged positive ions and 
the probability that the first cell contains 
ions and n,_ electrons, the second cell contains ions and electrons 
etc. the product the probability any specific combination the indivi- 
dual particles into the cells, which leads the required occupation numbers 
and the total number different combinations leading these same 
occupation numbers. 

Using these occupation numbers the probability density 
combination is, according equation (55), proportional 


(56) 


n m+n 
where and are the radius vectors the centers the and jth cells and 
where the contributions pairs particles situated the same cell are 
omitted from the right. This latter step justified the fact that the final 
results are independent the size the cells. 

The total number different combinations leading the required occupa- 
the total number electrons (or ions). Setting and 


for the required probability 


exq| - (6214 — | 
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From Stirling’s formula 


1 


then one may write and thus 


1 


since Substitution (58) into (57) results 
i 


Since the are very large numbers, one may regard and 
continuous variables and one may then consider equation (59) the probability 
density M24, ... space which must multiplied the volume 
element obtain the probability. Each side equation (59) 
can then also interpreted the probability that the number particles 
the first cell between and etc., since this leads volume 

The Jacobian this transformation constant that the pro- 


Another transformation introduced expanding and into 
Fourier series with periodic boundary conditions, assuming the total volume 
have the shape cube volume 
The equations transformation are then 


2 


since and are definition only determined points spaced apart 
the highest value each component taken the summations 
equation (61) then there are just enough equations determine the co- 


efficients and unambiguously. 
Since and are complex, convenient write 


and then (61) and (62) may rewritten 


where 
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and there are corresponding equations with the subscript instead 
convenient regard the transformation from space 
Space indicated (61) and (62) the sum the two transformations 
(63) and (64). The Jacobian (63) from space Sks 
space clearly constant and thus the volume element 
proportional the volume element The further trans- 
formation (64) can split into individual transformations which have the 
nature transformations from two-dimensional cartesian polar co-ordinates 
space then proportional 


(65) Tk A ka: Ay daygdays. 


remains rewrite equation (60) terms the Oks 
co-ordinates. convenient regard the for this purpose con- 
tinuous functions position defined equation (61) and rewrite equation 
(60) substituting integrals for the sums. 


The co-ordinate transformation relations (61) may written integrals 
space. the volume very small, the limits integration may taken 
infinity and the first equations (61) becomes the Fourier transform 


relation 


and similar equation applies Equation (67) may modified using 
the convolution theorem Fourier transforms, the form given equation 
(10) (11) Booker, Ratcliffe, and Shinn (1950), obtain 


;—T ;=const 


where kept constant the integrations indicated the left and the 
right. similar equation applies the Setting (68) results 


Similarly the first integral the curly bracket (66) may modified 
with the aid (68) and using the new set space variables 
instead the old set and 
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The inside integral (70) best evaluated changing spherical polar 


0 


This integral does not converge the limit taken infinity. Taking the 
limit temporarily where very large value |u,| one 
obtains 


kL 
0 


and thus (70) may rewritten the form 


The value the the right (71) does not change omitting the 
term coskL because the integral cos kLdk vanishes tends 
infinity. Equation (66) may then rewritten with the aid (69) and (71) 


Using equation (65) and (72) one can calculate the average value 


(73) 
and similarly 
(74) 2N/V. 
one defines analogy (61) 
one may write 
view the symmetry the distribution (55) and thus 
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should recalled this point that defined equation 
analogous (67) 
Writing for the electron density deviation the inverse form 
(77) may written 


Equation (76) shows that the volume taken unity then for sufficiently 
large values the mean square amplitude the Fourier component 
the electron (and ion) density fluctuations equal the electron density 
This value would exist random distribution electrons. the other 
hand then the mean square amplitude the Fourier com- 

Equation (73) shows that the volume taken unity, the mean square 
amplitude the Fourier component the difference between the electron 
density and the ion density for sufficiently large values This value 
the Fourier component again appropriate random distribution and 
general reduced the ratio Thus for the values 
for which where the Debye length, the Fourier 
component the electron density distribution almost equal the Fourier 
component the ion density distribution and the mean square value both 
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NOTES 


THE ELECTRICAL RESISTANCE MAXIMUM DILUTE MAGNESIUM ALLOYS 


recent publication Hedgcock al. (1960) some low temperature 
electrical resistivity measurements were reported dilute magnesium and 
aluminum alloys. the purpose the present note report the occurrence 
resistance maximum discovered alloys higher 
manganese content than the alloys described the earlier publication. 

All the resistance specimens were cut from 
ingot having nominal manganese content at.%. Resistance samples 
obtained from the ingot were quenched from various temperatures and the 
resistance ratio measured for each sample. Table contains values 
the resistance ratio and estimate the manganese content solid 
solution. The nominal concentration manganese given Table was 
deduced from the phase diagram the system (see Metals Handbook (1948)) 
and indicates that the resistance ratio increases with increasing manganese 
concentration would expected solid solutions were formed. order 
keep much manganese possible solution, the resistance specimens 
received strain-relieving anneal after being quenched and before being 
measured. 


TABLE 


Estimated atomic 


percentage 

0.152 0.16 1.2 21.3 17.8 
0.212 0.22 1.9 22.6 15.0 
0.328 0.35 3.3 25.4 10.0 
0.371 0.40 4.0 27.0 8.8 
0.498 0.56 6.3 31.0 5.8 

6.7 33.0 5.2 


0.527 0.60 


Figure shows some typical temperature resistance curves the region 
the resistance maximum. All the maxima were preceded pronounced 
minimum the electrical resistance. Table gives the temperature 
occurrence for both the resistance maximum and minimum well 
normalized value the the resistance minimum 
can seen that the position the resistance minimum moves 
higher temperatures the concentration manganese increased and that 
the normalized depth the minimum decreases with increasing impurity 
content. can also seen Table that the temperature the resistance 
maximum, increases the impurity content increased. 
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Mg-Mn 


Mg-Mn (~040 at . %) 
Mg-Mn WQI6 at. %) 


5 6 ' 2 3 4 5 6 8 9 
T (@K) T (@K) 


Fic. Temperature dependence the relative resistance number magnesium alloys 
described Table 


The anomalous resistance behavior the sample with nominal manganese 
content 0.16 at.% shows evidence levelling off even down 0.2° 
Preliminary measurements indicate negative magnetoresistance for these 
alloys the region the resistance maximum. Further details the electron 
transport and magnetic properties these alloys will given later and 
more detailed publications. 

The fact that manganese known induce resistance maximum, the 
monovalent metals gold, silver, and copper (see, for example, Korringa and 
Gerritsen (1953)) and the divalent metal magnesium* adds strength the 
argument that only the specific nature the impurity needs considered 
discussing low temperature resistive anomalies. 


The authors would like thank the Dow Chemical Metallurgical Labora- 
tories for supplying the magnesium alloys and the National Research Council 
Canada for the support this research. 
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*Magnetoresistance measurements made Muto (1959) have indicated the probable 
occurrence resistance maximum the divalent metal zinc. 
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LETTERS THE EDITOR 


Under this heading brief reports important discoveries physics may published. These 
reports should not exceed 600 words and, for any issue, should not later than six weeks 
previous the first day the month issue. proof will sent the authors. 


Frequency Measurement Standard Frequency 


Measurements are made Ottawa, Canada, using N.R.C. caesium-beam frequency resonator 
reference standard (with assumed frequency 770 c.p.s.). Frequency 
deviations from nominal are quoted parts per negative sign indicates that the 
frequency below normal. 


GBR, 


Date, 
May MSF, 6-hour 24-hour WWVB 
1960 60 kc/s average* average 60 kc/s 
—140 —142 N.M. 
—148 —147 N.M. 
N.M. N.M. N.M. 
—147 —144 —144 N.M. 
—124 —140 N.M. 
—142 —144 —144 
—143 —144 —142 —64 
—134 —137 —139 —65 
—146 132 —64 
N.M. N.M. N.M. 
N.M. N.M. N.M. 
—138 N.M. N.M. N.M. 
Midmonthly 
mean —141 —142 
Midmonthly 


mean of WWV — 137 


Note: N.M. no measurement. 
*Time of observations: 01.00 to 04.00 U.T.; 10.00 to 13.00 U.T. 


RECEIVED JUNE 17, 1960. S. N. KALRA 
DIvIsION OF APPLIED PHYSICS, 

NATIONAL RESEARCH COUNCIL, 

Ottawa, CANADA. 


l[ssued as N.R.C. No. 5746. 
2Cf. Kalra, S. N. 1959. Can. J. Phys. 37, 1328. 
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MANUSCRIPTS 


General.— Manuscripts, English French, should typewritten, double spaced, 
paper in. The original and one copy are submitted. Tables and captions 
for the figures should placed the end the manuscript. Every sheet the manuscript 
should numbered. Style, arrangement, spelling, and abbreviations should conform the 
usage recent numbers this journal. Greek letters unusual signs should written plainly 
explained marginal notes. Characters set boldface type should indicated 
wavy line below each character. Superscripts and subscripts must legible and carefully 
placed. Manuscripts and illustrations should carefully checked before they are submitted. 
Authors will charged for unnecessary deviations from the usual format and for changes 
made the proof that are considered excessive unnecessary. 

Abstract.—An abstract not more than about 200 words, indicating the scope the work 
and the principal findings, required, except Notes. 

References.—References should listed alphabetically authors’ names, un- 
numbered, and typed after the text. The form the citations should that used current 
issues this journal; references papers periodicals, titles should not given and 
only initial page numbers are required. The names periodicals should abbreviated 
the form given the most recent List Periodicals Abstracted Chemical Abstracts. All 
citations should checked with the original articles and each one referred the text 
the authors’ names and the year. 

Tables.—Tables should numbered roman numerals and each table referred the 
text. Titles should always given but should brief; column headings should brief and 
descriptive matter the tables confined minimum. Vertical rules should not used. 
Numerous small tables should avoided. 


ILLUSTRATIONS 


figures (including each figure the plates) should numbered con- 
secutively from up, arabic numerals, and each figure referred the text. The author’s 
name, title the paper, and figure number should written the lower left corner the 
sheets which the illustrations appear. Captions should not written the illustrations. 

Line should carefully made with India ink white drawing 
paper, blue tracing linen, co-ordinate paper ruled blue only; any co-ordinate lines that 
are appear the reproduction should ruled black ink. Paper ruled green, yellow, 
red should not used. All lines must sufficient thickness reproduce well. Decimal 
points, periods, and stippled dots must solid black circles large enough reduced 
necessary. Letters and numerals should neatly made, preferably with stencil (do NOT use 
pewriting) and such size that the smallest lettering will not less than high 
when the figure reduced suitable size. Many drawings are made too large; originals 
should not more than times the size the desired reproduction. Whenever possible 
two more drawings should grouped reduce the number cuts required. such groups 
drawings, large drawings, full use the space available should made; the ratio 
height width should conform that journal page in.), but allowance must 
made for the captions. The original drawings and one set clear copies (e.g. small 
photographs) are submitted. 

Photographs.—Prints should made glossy paper, with strong contrasts. They 
should trimmed that essential features only are shown and mounted carefully, with 
rubber cement, white cardboard, with space between those arranged groups. mount- 
ing, full use the space available should made. Photographs are submitted 
duplicate; they are reproduced groups one set should mounted, the duplicate 
set unmounted. 
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